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Abstract 

We investigate the spectacular collider signatures of macroscopically displaced, neutral par- 
ticles that decay to Higgs bosons and missing energy. We show that such long-lived particles 
arise naturally in a very minimal extension of the Standard Model with only two new fermions 
with electroweak interactions. The lifetime of the long-lived neutral particles can range from 
10 -2 mm to 10 6 mm. In some regions of the parameter space, the exotic signals would have 
already been selected by the ATLAS and CMS triggers in their 7 and 8 TeV runs, hence hiding 
in the existing data. We also discuss the possibility of explaining the mild anomalies observed 
in the diphoton Higgs channel and the WW production at the LHC. 
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1 Introduction 



The recent discovery of a new resonance with mass ^125 GeV [1] at the Large Hadron Collider 
(LHC) is an extraordinary milestone for fundamental physics. Although the properties of the 
resonance are consistent with those of the Higgs boson of the Standard Model (SM) at the 2a level, 
the experimental uncertainties are still large and allow ^ 10 % deviations from the SM, even if we 
put aside the potential indications of a mild excess in the diphoton channel, pp —> h —> 77 [2]. If 
the diphoton excess or any other deviation in the Higgs properties is indeed confirmed, it means 
the presence of new physics beyond the SM at the electroweak scale. 

In this paper, we explore the possibility of new physics right at the 100- GeV scale providing 
highly non-standard Higgs signals that elude existing search strategies aimed at the SM Higgs or 
its close variants. Specifically, the Higgs boson can be produced in association with missing energy 
from the decay of a new, macro scopically long-lived, neutral particle produced via electroweak gauge 
interaction. In some regions of the parameter space, it is even possible that the signals we propose 
are already in the existing LHC data, while in other regions the near-future discovery potential is 
high. Our work differs from most of the existing literature on non-standard Higgs phenomenology, 
which have either focussed on exotic Higgs decays [4, 6, 7, 8], or prompt Higgs production from 
new physics, for example, in the context of Gauge-Mediated Supersymmetry Breaking (GMSB) 
[11, 12], though displaced Higgs scenario has been briefly discussed in [9, 10] 

The extension of the SM proposed in this paper contains only two new fields: a massless 
singlet fermion and a massive electroweak-triplet fermion. 1 We do not introduce any new gauge 
interactions beyond the SM. Our model thus significantly differs from the hidden valley models [3], 
which involve new confining gauge interactions, although a class of hidden valley models can also 
give rise to non-standard Higgs signals similar to the ones we propose [4]. The lagrangian of our 
model and its field theoretic properties will be discussed in detail in Section 2. 

Despite its minimality, our model displays a strikingly rich collider phenomenology, with a 
two-dimensional parameter space divided into eight signal regions, depending on the spectrum and 
dominant decay modes. In Section 3, we will define these signal regions and analyze the constraints 
and discovery potential for each region. We discuss how our model could lead to a diphoton excess, 
and also explain the mild anomaly in the sample observed at the LHC [31] in a manner 

similar to the charginos in supersymmetric models [24]. 

2 The Model 

We consider a model with additional spin- 1/2 fermions that are vectorlike under the SM gauge 
group. We introduce a Dirac fermion uj and a left-handed Weyl fermion Xol with SM gauge quantum 
numbers (1, 3)o and (1, l)o, respectively. We assume a new unbroken global U(l) symmetry, U(l) x , 
under which uj and Xol both carry charge +1, 2 the importance of which will be explained shortly. 

1 These new states could further be embedded in a supersymmetric model if one wishes, although we will not go 
through such model building exercise in this paper. 

2 U(1) X is anomalous in the presence of gravity, but the anomaly could easily be cancelled by simply adding extra 
SM singlet fermions with U(l) x charge, as there are no anomalies induced by SM gauge interactions. 
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Figure 1: A schematic diagram showing the mass hierarchies of the new fermions in our model with 
rriu = 140 GeV for the two cases: (i) A 2 < 200 TeV, (ii) A 2 > 200 TeV. The exact mass-splittings 
primarily depend on the scale A 2 and to a much lesser degree on the scale A3 as long as A3 > lOTeV. 

After electroweak symmetry breaking (EWSB), we will use the electric-charge basis defined as 

cj 1 — icj 2 . id 1 + icj 2 o 

= — 7 =- ) U _ = —j=-, U = U> , (1) 

where uo a (a = 1, 2, 3) label the components of the SU(2) L -triplet uo. The three fermion fields 
uo f _ and uo carry electric charge +1, —1 and 0, respectively, while they all have U(l) x charge +1. 
The antiparticle of UJ+ will be referred to as cj_, which is distinct from uo'_, as the former carries 
U(l) x charge —1 while the latter +1. As we will see below, uj± and gj' t may even have different 
masses. The antiparticle of uo (carrying U(l) x charge —1) will be referred to as cJo- Similarly, the 
particle and antiparticle interpolated by the Xol field will be denoted by xo and Xch respectively. 
The most general renormalizable lagrangian £ ren consistent with SM gauge invariance and the 
U(l) x symmetry is given by 

£ ren iuJ^B^u -m^uJuj^r ixl^d^X^ • (2) 

Here, the new fermions have only gauge interactions. In particular, they have no direct couplings 
to the Higgs or SM fermions, and Xol is completely decoupled. The U(l) x symmetry is crucial 
here; without it, we would have a dangerous renormalizable interaction H^(j a cr 2 l\^uj^ with the SM 
Higgs field H and lepton doublet ^ L , which would mix the new fermions and SM leptons. To avoid 
severe constraints from lepton flavor violation, we have chosen uo to be Dirac rather than Majorana, 
with the conserved charge, U(l) x . 

The absence of direct couplings between the new fermions and the Higgs is an artifact of 
renormalizability, and hence vulnerable to nonrenormalizable interactions. At dimension-5, three 
non-gauge interactions arise between the new fermions and the Higgs field: 

0W=H*H<4!cjZ, C?f=ie a6c ffW w 6t w c j of =rfa a H^X^, (3) 
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where a a are the Pauli matrices. Again, the U(l) x symmetry plays a crucial role here in forbidding 
dangerous operators like H T a 2 a a H e^u®. The absence of direct coupling between the new and SM 
fermions therefore persists at dimension-5. There are three more dimension-5 operators involving 
the new fermions and SM gauge field strengths: 

{ 5 4) = g 1 B IM/ c^VaX . Of = ig 2 e abc W« v co b Ja^co c L , of ] = g 2 W^ c^Wx^. . (4) 

The full Lagrangian we consider is given by 

Q(n) 

C = C SM + £ren + V -f- (5) 
n An 

with A n > 1 TeV. 

(2) (3) 

We focus on the effects of 0\ J and 0\ J for all of our studies of the LHC phenomenology unless 
otherwise noted. O^p just gives a universal mass shift of order v 2 / Ai to upon EWSB, and 
induces tiny couplings, suppressed by Ai, between uj and the Higgs boson /i, the signals from which 
would be swamped at the LHC by those from the renormalizable gauge interactions unless Ai is as 
low as the TeV scale. We will therefore ignore hereafter unless otherwise noted. In contrast, 

(2) (3) 

the operators O b and 0§ upon EWSB induce mass splittings between different components of 
the triplet cj, thereby opening up phase space for transitions between those states via a virtual W ± . 
Specifically, the operator 0\ ) after EWSB (v ~ 246 GeV) gives mass terms 

v 2 + + 

thereby splitting the masses of the two charged states uj± and cj^ by ~ v 2 where the splitting 
can be as large as several GeV for A2 > 10 TeV. The neutral component cjq receives no contributions 
from this operator. Like those from the Yukawa couplings from between uj and h would 
be completely swamped by the renormalizable gauge interactions, and hence not interesting. On 

(3) 

the other hand, the operator 0\ J upon EWSB induces a mass term 

v 2 + 

"2A^^0L + h -c, (7) 

thereby slightly mixing the neutral states c^ol and Xol 5 while leaving the charged components 
unchanged. The mixing induces a Yukawa coupling of the size v/As among cjor, Xo and /i, thereby 
dramatically altering the phenomenology of xo 5 as it was completely inert at the renormalizable 
level. The increase of the cjq mass due to this mixing is of order v 4 / 'A^rn^,, which can be at most 
several tens of MeV for A3 > 10 TeV and ^100 GeV, and hence its phenomenological effects are 
quite minimal. Therefore, we expect a rich collider phenomenology in a two dimensional parameter 
space spanned by A2 and A3, where the former controls the charged sector by creating a mass 
gap between uo± and cj^, while the latter governs the neutral sector by giving a uo-xo-h Yukawa 
coupling. In particular, the smallness of these mass splitting and Yukawa coupling gives rise to the 
interesting possibility that the transitions between the new fermion states involving h and/or W ± 
can occur with possibly macroscopic decay lengths. 
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We will ignore the dipole operators C^ 4 ' 5 ' 6 ^, with the following justifications. First, the following 
very simple UV completion of the effective Lagrangian (5) suggests that it is natural to take 
^4,5,6 ^> Ai,2,3- Imagine a heavy Dirac fermions ^ with SM gauge quantum numbers (1, 2) 1 / 2 and 
a U(l) x charge +1. The most general renormalizable lagrangian involving ip is given by 

i^D^ - - (A .hVc^Vl + A' ^ R u*a a H + A" ^XqlH + h.c.) , (8) 

where the SU(2) L -doublet indices among a a , ij) and H are implicit. Upon integrating out the heavy 
fermion -0, O^ 1 ' 2 ' 3 ^ are generated at tree level with Ai = A2 = M/AA' and A3 = M/AA". On the 
other hand, the dipole operators C^ 4 ' 5 ' 6 ^ are generated at 1-loop with A^^ ~ 167r 2 Ai 5 2,3, respec- 
tively. Thus, (9 5 ' ' are suppressed by two more orders of magnitude than (9 5 ' ' in the effective 
Lagrangian (5). The effects of C^ 4 ' 5 ^ are then completely swamped by those of the renormalizable 
gauge interactions, so we will ignore C^ 4 ' 5 ^ hereafter. The operator could induce a potentially 
interesting process uo —> xo + ~f/Z, 3 although with the (167T 2 ) 2 suppression, it would be completely 
subdominant relative to the 0\ -induced process, ojo — > xo + unless we (almost) close the phase 
space for the latter by having m w nearly at or below m^. Since we are interested in robust on-shell 
Higgs production in this paper and will take m w sufficiently above m^, we will also ignore 
hereafter. 

We are concerned with small mass splittings; therefore it is important to take into account 
mass splitting effects induced by and Z loops, which arise even at the renormalizable level. 
In the absence of nonrenormalizable operators, cj±, u f T and uo are degenerate at tree level by the 
SU(2) L symmetry. However, since electroweak symmetry is broken, we expect that the degeneracy 
should be lifted by loops. The one-loop radiative mass splittings between cjo, u± and from the 
renormalizable gauge interactions alone are given by [15, 16]: 



Am w = m u + - m^ n m,„/ - m 

OL2 



4vr^ 



where 

IT 



f{r) = r 4 log r + r (2 + r 2 ) y 4 — r 2 Arctan I 



(9) 



(10) 



Here, Am w is always positive, rendering ojo lighter than uj± and uj'^ (in the absence of the contri- 

(2) 

butions from 0\ } ). Recently, Am w has been calculated at two- loop level [17], which significantly 
reduces the renormalization scale dependence compared to one-loop computation. Numerically, the 
mass splitting varies between ~ 150-165 MeV for between ~ 100-1000 GeV. 
To summarize, the phenomenological Lagrangian of our model is 

£> — £>SM + Ami + ^mass + £gauge + 5 (H) 



3 Prompt and displaced j/Z signatures also arise in GMSB model [11, 13, 14]. 
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where, at the leading non-trivial orders in v/Az,3, 

£ kin = iu +1 »d^ + + iu;'_^d^'_ + icJo^d^uo + ixJ L ^^X 0L , (12) 

£mass = -m UJ± UJ+LO+ - w'_uj'_ - m Uo UJ UJ , (13) 

Asauge = 92 [W+ (w + 7^o - ujq^uj'J) + h.c] - g 2 cos(9 w Z^(uj + Yu+ - uj'^ujL) 

- eA„ (uj +1 »uj + - uj'_j"u'_) , (14) 



v 2 



with 



£a = ~^^r K ~ + h.c. - -^«X 0L + XJ L ^ 0R ) , (15) 



V 



it v _i_ vt 



A; 
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v 2 . i; 2 i; 4 



^± = ^ + Am w - — — , m,/ =m w + Am w + — — , m UQ = + ■ (16) 

Here, we have chosen A3 to be real and positive without loss of generality by adjusting the phase 
of Xol- The phase of A2 cannot be removed by field redefinition, however, and it would lead to 
CP- violating Yukawa couplings between uj and h. But since it would be hard to observe such CP 
violation effects, we have taken A2 to be also real and positive in Eq. (16) for simplicity. 

Before we move on to the discussion of collider phenomenology of the Lagrangian (11), we would 
like to make a few comments on possible variations of the model. One simple variation is to get 
rid of U(l) x and make uj Majorana instead of Dirac, though this possibility has flavor issues as 
discussed earlier. Another possibility is to keep U(l) x and make Xo a l so Dirac by adding Xor 1° 
the theory. These two variations have difficulties in cosmology, however. In the Majorana case, 
fine-tuning would be necessary to make xo massless, while in the Dirac xo case, there is no reason 
for xo to be massless. A massive xo would lead to over-closure of the universe unless its mass is 
< 1 eV, since xo is extremely weakly interacting and stable. In fact, the mass must be <C 1 eV, 
because if it is near ~ 1 eV, (a significant fraction of) dark matter would be hot, thus at odds with 
observation. Therefore, the simplest approach is to prevent xo from acquiring mass by symmetry 
as in our model. 

The big-bang nucleosynthesis (BBN) bound on extra relativistic degrees of freedom do not 
necessarily constrain the massless xo? because the xo gas is not in equilibrium with the SM gas 
unless the temperature is far above the weak scale. Thus, the temperature of xo can be much lower 
than that of the SM during BBN. Also, between the weak scale and the BBN temperature, the 
SM gas gets heated up several times, possibly through electroweak phase transition, and certainly 
via the annihilations of massive particles such as t, Z/W ± 1 6, etc., as the temperature drops below 
their masses. The xo g as > on the other hand, remains unheated. 



3 Collider Phenomenology 

In this section, we discuss the collider phenomenology of the lagrangian (11). We first study the 
production mechanisms of the new fermions in our model and their decay channels and widths. 
In particular, we emphasize that long-lived particles decaying to Higgs plus missing energy is a 
prominent signature of our model in a large portion of the parameter space. We will then divide 
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our parameter space into eight regions depending on the signatures, and discuss for each region 
the current experimental constraints on our model as well as discovery prospects. We will see that 
some of the exotic events would have already been selected by the ATLAS and CMS triggers in 
their 7 and 8 TeV runs and can be used for possible discovery. 

3.1 uo production, decay, and the benchmark point 

The new fermions uj can be Drell-Yan produced at hadron colliders through gauge interaction in 

£gauge E( l- ( 14 ) aS: 

PP ~ > l/Z* — * ^+ UJ - ? PP ~~ ^ l/Z* ~~ ^ u>_ , 

pp —> W + * — > uj+ cJo , PP — >• ^+ , 

pp — >> VF~* — y ujq uj— , pp — » W~* —> ujLujq • 

(At e + e~ colliders, only the 7/^* mediated processes are possible.) Recall the U(l) x symmetry 
under which cj + , uj'_ and ujo carry charge +1, while their antiparticles uj f + and uJo carry — 1, 
thus forbidding the final states like ujouj+ as the SM particles are all U(l) x neutral. Although, 
direct xo production at the colliders is highly suppressed by v/A%, the fermions uj eventually decay 
to xo, which is the lightest stable particle (LSP) in our model and appears as missing energy in 
the detector. 

Clearly, the signatures of our model depend on the decay of uj fermions. In Eq. (16), we observe 
that the effect of the scale A3 on the mass of the neutral particle ujo is at most a few tens of MeV 
for A3 ^> 10 TeV and ~ 100 GeV. Therefore, the mass hierarchy is more or less determined by 
the scale A2, with uj'^ being the heaviest among the three uj fermions 4 while the next-to-lightest 
stable particle (NLSP) could be either ujq or uj±. A schematic diagram for the mass hierarchies of 
the uj fermions is shown in Fig. 1 for = 140 GeV, which will be our benchmark point for the 
rest of the discussion unless otherwise stated. 5 Led by this observation of mass hierarchy among 
the uj fermions, we can roughly divide the parameter space of interest into two categories: 

(I) uj± NLSP : 

For A2 < 200TeV, the mass formulas (16) tell us that the V/A2 effects, which tries to lower m UJ± 
below m Wo , is greater than the radiative mass splitting Am w , which tries to push rn u± above 
777^ . Consequently, uj± is the NLSP with mass-splitting m Wo — m UJ± as big as a few GeV for 
A2 < 20 TeV. The only decay channel for the NLSP is uj+ — )► xo W + (or the charge-conjugate 
process uj- —> Xo W~) with decay width 

r K -> xo w+) = ™„ ± (17) 

On the other hand, the neutral uj has two decay modes available: 

4 For A3 < 10 TeV, it is possible for the charged fermions, u± and lj± to be lighter than the neutral cjo- However, 
in this case, all the uj fermions decay promptly to xo, since is not small. This signature is phenomenologically 

no different from the Region (iv) signals to be discussed in section 3.3. 

5 A reason for this benchmark point is to have ra^Q > mh such that the Higgs boson h in ujo — » xo h is always 
produced on-shell. 
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(i) loq — > xo h (or the charge-conjugate process ujq — > Xo h) with a decay width 

,2 \2 



iP 1 { 771? \ 



(18) 



(ii) cjo — > ^+ * (or the charge- conjugate process ujo — >> uj- with 

2 / m 2 \ 

r(cj cj + Pi) = ^3 g f ^ ^ [j^2 J Q ( 6m " - m i) > 

T(uj° — )► cj + + hadrons) 

'o ;k w <m,+ (19) 

< r^o^ +7r - 5m w ^1 - J ; ^tt+ < ^o; < 

where Sm u = m^ — ra^ ± , iV c = 3, ~ 92.2 MeV is the pion decay constant, Gp is the 
Fermi constant, and the function F is given by 



F{a) = \/l - a ^1 - - 4a 2 ^j + y a 2 Arctanh^/T 3 " 



(20) 



The coefficient ft appearing in Eq. (19) is an 0(1) factor introduced to parametrize the 
uncertainty in the hadronic decay width of cj, which we estimate to be on the order of the 
color factor N c times the leptonic decay rate, r^^^ e - i?e . In our analyses below, k will 
be varied between and 2. 

The heaviest of the new fermions uj' t can decay to lighter states cjo, cJo 5 Xo and Xo through 
the channels, uj'_ — >> xo W~ , uj'_ —> ujo W~* and their charge- conjugate processes. Their decay 
widths can be calculated by the replacements m LJ ± in Eq. (17), and ujo —> uj'_ and 

uj+ — >► cjq in Eq. (19). Note that the direct decay of uj'± to uj± is forbidden by U(l) x symmetry 
(but it can happen in two-step weak decays with an intermediate ujq). 

(II) uj NLSP : 

If the scale A 2 > 200 TeV, the mass splittings between the uj fermions are dominated by Aui^. 
In this scenario, the NLSP is ujo, while uj± and uj' t are nearly degenerate with each other 
for A2 ^> 10 3 TeV. The only decay channel of the NLSP is ujq — » xo h, while the fermions uo± 
and uj' t have additional decay modes uj+ — » ujo W + and uj'_ ujq W~ besides uj+ xo W + 
and uj'_ — >> Xo^~5 respectively. The decay widths for cjo, cj± and cj^ can be deduced from 
Eqs. (17)-(19) by appropriate replacements as in the previous case. Since uj' t is nearly degen- 
erate with cj±, its decay widths are the same as the corresponding widths of uo±. 

In Fig. 2, the branching ratio and decay lifetimes of the uj fermions are plotted in the A2-A3 
parameter space for our benchmark point. A few comments are in order: 

• The lifetime of the heavier charged fermion a/p shown in Fig. 2(c), never exceeds « 100 mm. 
Naively, one might expect the lifetime to increase as the scale A3 increases and suppresses 
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Figure 2: Top Row: Branching Ratio (BR) contour curves are plotted for (a) cjo, (b) uj± and (c) u f T 
fermions in the A2-A3 parameter space for the benchmark point = 140 GeV. The central solid 
red line depicts a BR = 50 % for the two decay channels. The two solid black lines correspond to 
BR = 99 % for each of the two decay channels. Finally, the dashed lines for each BR contour line 
depicts the uncertainty in the calculation of hadronic decays by varying the n parameter between 
0-2 (see text for definition of n). Bottom Row . Lifetimes of the uj fermions as contours in the A2-A3 
parameter space for = 140 GeV. 



the W ± -xo-(^ , - = f coupling Eq. (15). However, since the splitting between ojq and oj^ remains 
> 150 MeV due to the radiative contribution Eq. (9), there is always a phase space for the 
weak decay uo f _ —¥ uo VF - *, keeping the u f T lifetime below « 100 mm. The leptons/hadrons 
from this W~* will be too soft to be detected at the LHC or Tevatron. Consequently, the uj' t 
production will just appear as the uo production for all practical purposes. 

• The same, however, does not hold true for the lighter charged fermion uo±. Towards the 
upper-left corner of Fig. 2(b), uj± becomes the NLSP, with the only possible decay mode being 
cj + —> xo W + . The lifetime of this decay increases as A3 increases and suppresses the cj±-xo- 
W T coupling Eq. (15), letting the uo± travel deep into the detector with lifetime > 500 mm. 
Towards the upper-right corner of Fig. 2(b), the NLSP becomes ojo with the mass splitting 
between ojq and uo± dominated by the radiative effect Eq. (9), i.e., Tri(jj_^_ 772 r^> 150 MeV. 
Here, like the uiL case above, the u>± production will be equivalent to the coq production for 
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all practical purposes. 

• The neutral fermion, ujq can be comparatively long-lived as shown in Fig. 2(a). This follows 
from the fact that ujo is the NLSP in majority of the parameter space (A2 > 200 TeV) and 
hence the only decay channel available is ujq — >> xo h. For ujq lifetimes exceeding ~ 10 mm, 
its decay results in Higgs production amidst the detector leading to a plethora of collider 
signatures as will be discussed next. 

3.2 Simulation 

The events for uj pair production were generated in MadGraph [18] and showered in Pythia [19]. 
The input model file for MadGraph was generated using the Mathematica package FeynRules [20]. 
The uj pair-production cross-sections, shown in Table 1, were obtained using MadGraph. 





>/s = 1.96 TeV 


v^ = 7TeV 


= 8 TeV 


y/s = 14 TeV 




Tevatron Run 


LHC Run 


LHC Run 


LHC Run 


pp(p) -> U UJ- 


0.16 pb 


0.61 pb 


0.78 pb 


2.01 pb 


pp(p) -> 


0.16 pb 


1.23 pb 


1.52 pb 


3.45 pb 


pp(p) -> u uj' + 


0.16 pb 


1.23 pb 


1.52 pb 


3.45 pb 


pp(p) -> LOqJ_ 


0.16 pb 


0.61 pb 


0.78 pb 


2.01 pb 


pp(p) — > W+ 0J- 


0.22 pb 


0.91 pb 


1.14 pb 


2.74 pb 


pp(p) — > Uj' + Uj'_ 


0.22 pb 


0.91 pb 


1.14 pb 


2.74 pb 



Table 1: Cross-sections for pair-production of uj fermions at the LHC and the Tevatron for the 
benchmark point m u ~ 140 GeV. 

While the proper lifetimes of the uj fermions has been plotted in Fig. 2, the actual decay length 
inside the detector will, event by event, be different from the proper lifetime depending on the boost 
factor. To study this effect, we have performed an elementary detecter simulation to calculate the 
fraction of uj particles that decay in different parts of the LHC detector. The lifetime of a particle 
in the lab frame, r is related to its proper lifetime, tq by 

^71^™- (21) 

The mean transverse distance Lx travelled by the particle in the lab frame is 

L T v T r v T -fr —T . (22) 
m 

The pre-factor in the above equation p T /m is precisely the distribution plotted in Fig. 3(a). If TVo 
is the number of long-lived particles produced at the primary vertex, the fraction N(x T ) of these 
particles decaying within a distance x T in the transverse direction is given by 

N(x T ) = (l- e~ x T /LT )N (23) 
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8 TeV LHC , A 2 = 10 3 TeV 




Figure 3: (a) p^/m distribution of cjq for 8TeV LHC run for = 140 GeV and A 3 = 10 3 TeV (i.e., 
the cjq NLSP region). The normalization on the vertical axis is arbitrary, (b) Fraction of cjq decaying 
inside the various components of the ATLAS detector as a function of A3 for m w = 140 GeV and 
A2 = 10 3 TeV. Only decays with transverse decay length > 0.1 mm are included in the plot. 



The fraction of cjq decaying inside different sections of the ATLAS detector are shown in Fig. 3(b) 
for A3 » A2 (i.e., the uo NLSP region). Although our simulation has used the dimensions of ATLAS 
detector, the results also apply to the CMS detector to a good approximation. Note that uo with 
lifetimes as long as ~ 10 m can have significant decays inside the inner detector. 

3.3 Constraints and LHC Discovery Prospects 

This section will discuss the current experimental constraints on our model as well as future discov- 
ery prospects. For clarity, we classify the parameter space of our model into eight different regions 
based on the final states as shown in Fig. 4: 

Region (i): Long-lived charged particles 

Here, the charged fermion uj± has lifetimes exceeding 50 cm. Current LHC searches for long-lived 
massive charged particles are sensitive to this scenario, triggered by activity in the muon tracker 
from long-lived charged particles or by the large missing energy utilizing only the calorimetric activ- 
ity. The most constraining of these long-lived charged particle searches is by the LHC experiments 
[32] at = 7 TeV. In particular, for lifetimes long enough for the uj± to reach the muon system 
of the detector, the parameter space is tightly constrained. For example, in the ATLAS search, the 
signal efficiency times acceptance for a O(100 GeV) long-lived charged particles that pass through 
the detector is ~ 20%. From our elementary detector simulation shown in Fig. 5 for two represen- 
tative points in the region of parameter space of interest, we estimate that the benchmark point of 
= 140 GeV yields > 0(1) events for A 3 > 2 x 10 7 TeV. Another ATLAS search [33] has looked 
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A 2 [TeV] 



( i) Long-lived charged particles 

(ii) Displaced WW + MET 
(Hi) Prompt WW + MET 

( iv) Prompt W h + MET and 
Prompt WW + MET 

(v) Displaced W h + MET and 
Prompt WW + MET 

( vi) Displaced W h + MET and 
Displaced WW + MET 

( vii) Displaced hh+ MET 
iyiii) Pure MET 



Figure 4: The A2-A3 parameter space classified according to the most dominant final states of the 
uj decays at the LHC for m w = 140 GeV. The boundary between two given regions is defined as 
points where the two corresponding final states become comparable to each other. 



for long-lived charginos with disappearing track signatures. This search is sensitive to charged 
particles that decay in the outermost parts of the inner detector, i.e., with decay lengths ~ 50 
cm. Combined with other long-lived charged particle searches [34], we conservatively conclude that 
our benchmark point = 140 GeV is excluded in Region (i), although a more thorough detector 
simulation is warranted. 

Region (ii): Displaced WW + $^ (diagrams in Fig. 6) 

This region of the parameter space is characterized by the following two features: 

• uj± is the NLSP. uo and u f T promptly decay to uj± through weak interactions with the decay 
products from W^ being too soft to be detected. 

• cj± is long-lived, but unlike Region (i), its decay length is in the range between ~ 1 mm and 
~ 50 cm. Below 50 cm, direct detection of uj± at the LHC rapidly becomes challenging due 
to an insufficient number of hits in the inner detector. 

Therefore, the only potentially observable final state is the + from the long-lived uj± 

decays, where both W ± are hard and displaced from the primary vertex. The hadronic decays of 
the displaced W ± are difficult to observe due to large QCD background. For the leptonic decays of 
the W ± 1 while triggers for £ + |? T already exist at the LHC, we cannot reconstruct the secondary 
vertex from which the lepton originates. However, searches for kinked tracks could be a possible 
mode for discovery in this region of parameter space. 

Region (Hi): Prompt WW + |? T (diagrams in Fig. 6) 

As in Region (ii), uj± is the NLSP here, and the only relevant final state is VF + VF~ + with 
other decay products from W^ being too soft to be observed. Unlike Region (ii), however, uj± 
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Figure 5: Fraction of uj± decaying inside the various components of the ATLAS detector for 
rriu = 140 GeV for (a) A 2 = 100 TeV and (b) A 2 = 10 3 TeV. Only decays with transverse de- 
cay length xt > 0.1 mm are included in the plot. 



decays promptly, with lifetime < 1 mm. Ignoring the soft decay products, we observe that the 
Feynman diagrams shown in Fig. 6 are topologically similar to chargino/neutralino pair production 
in supersymmetric models. We have checked that our model is not excluded by opposite-sign 
dilepton searches [28, 29] since they usually require hadronic energy or much larger missing energy 
compared to our signal. The rriT 2 searches in [28] considered charginos decaying to sleptons which, 
in turn, give the dilepton signature. It is not immediately clear if this search is applicable to 
our scenario, though an optimized my 2 search could hold discovery prospects in this region. The 
trilepton searches [27] are moot for our scenario, since we have at most only two hard leptons. 
Finally, notice that the U(l) x symmetry implies that the final, on-shell W pair must come in 
opposite-sign (see Fig. 6). Therefore, this region of parameter space is also safe from same-sign 
dilepton searches [29, 30]. 

The + $t production in this region will also manifest as an excess in the SM WW 

searches [31]. Such an excess can actually improve the agreement with the experiment compared 
to SM contribution alone, as demonstrated in [24] for the case of WW production from ~ 110 GeV 
charginos decaying to Ws and gravitinos. The analysis in [24] also shows that, having similar 
kinematic features to the SM WW production, our + |? T are unlikely to contaminate the 

h — >► WW searches, unlike the new physics of the type discussed in [23], for example. Therefore, 
while the uncertainties in the current measurement are too large to make definitive statements, the 
2£ + $t searches could be an interesting path to discovery should our model prove to be correct. 

The operator which is mostly inconsequential and has thus been ignored so far, can have 

an interesting effect in Region (Hi). Although outside the plot in Fig. 4, the Higgs diphoton rate, 
pp (h —> 77) + X, can be significantly modified if Ai ~ A 2 ~ 1 TeV. Including the contributions 
of the top quark, W ± and the new charged fermions, ou± and a/p, the decay width of the Higgs to 
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Figure 6: Feynman diagrams leading to + j5 T final states in Regions (ii) and (Hi) of the 

parameter space, where uj± is the NLSP. Particles too soft to be detected at the LHC are shown 
in orange. Observe that, in each diagram, the on-shell VF-pair have opposite signs. 



diphoton is given by 



T(h -> 77) 



128\/2 7r 3 



A\{t w ) + -A h l/2 ( Tt ) + 9 hf+f _-^A h 1/2 (T f± ) 



f=Ul,L0' 



(24) 



where ti — m|/(4m?), the functions Ay 21 A\ are defined in [49], and the couplings ghf+f- f° r the 
charged uj fermions induced by the operators C^ 1,2 ^ are given by 



9h 



V V 

Ai A 2 



v v 
Ai A 2 



(25) 



For simplicity, we consider the UV completion (8) so that Ai = A 2 = A, where A in general is 
complex. The diphoton rate i? 77 is plotted in Fig. 7, where 



a(pp -> h + X) x Sr(/i ->> 77) 



a(pp^ h + X) x Sr(/i^77)| SM ' 
We find that the diphoton rates are maximized if A is close to real and positive. 



(26) 
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Figure 7: Contour plots of the diphoton rates R 77 (dashed black lines, white boxes) and the masses 
of lightest uj fermion, m^ (solid red lines, yellow boxes) are shown for A3 ^> Ai,A2, assuming 
Ai = A2. In plot (a), contour lines are plotted in Ai-ra^ plane for real and positive Ai } 2- In 
plot (b), the magnitude of the scales Ai and A2 is fixed at 1 TeV while their complex phase Arg[A] 
is varied. 



Region (iv): Prompt Wh + ^ T , Prompt WW + $ T (diagrams in Fig. 8) 

In this region, all the uj fermions decay promptly to xo in a single step. Fig. 8 shows two possible 
final states of interest at the LHC, Wh + $ T and W + W~ + $ T . Undoubtedly, the strongest 
constraint for this final state comes from the associated Higgs production searches in the 66 channel, 
pp -> W + (h -> bb). The 95% CL limits on a(pp -> Vh) x BR{h -> bb) (V = W, Z) by ATLAS 
and CMS for ~ 125 GeV Higgs are 1.8 and 2.5 times the SM value, respectively [26]. Most of our 
signal events pass the kinematic cuts used in these searches for < 200 GeV, thus comfortably 
excluding the benchmark point m w = 140 GeV. This situation is similar to a supersymmetric 
scenario with fi ^> Mi,M2, where the NLSP neutralino can decay to Higgs and the LSP, and the 
chargino decays to W and the LSP [21, 22]. 

Region (v): Displaced Wh + ^ T , Prompt WW + $ T (relevant diagrams in Fig. 8) 

In this region, the neutral fermion cjq is long-lived with lifetime in the range 1 mm-1 cm. The 
charged states, uj± and cj^, still decay promptly with < 1mm. This difference stems from the 
proximity of rrih to m w , which leads to a slightly more suppressed phase space for uo —> hxo than 
uj-\. —> W^xoi while the interaction strengths behind these decays are parametrically the same. 
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Since the lifetime of ujq is < 1 cm, it typically decays to h + ^ T before reaching the inner detec- 
tor, so the actual final states to be observed are the B hadrons from the h. 6 Notice, however, that 
these B hadrons originate from a vertex that itself is already displaced, as the decay ujo — >> h + ^ T 
has lifetime > 1 mm. This significantly degrades the efficiency of standard 6-tagging algorithms [48] 
that use a "signed impact parameter" as a discriminating variable with the positive sign being pre- 
ferred. As illustrated in Fig. 9, the sign is determined by the angle between the decay length vector 
and the jet-axis, where the sign is taken to be positive if this angle is less than 90°. As shown in 
Fig. 9(b), a significant fraction of b jets originating from displaced ujo decays will give negatively 
signed impact parameters, resulting in a much reduced efficiency of tagging the b quarks in our 
signals, although a detailed simulation of this effect is beyond the scope of a theoretical paper. Let 
us conclude that, unlike Region (iv), we expect that the Wh searches should not be very constrain- 
ing, and, turning this around, the observation of events with negative impact parameters should 
be regarded as the opportunity to probe this scenario. 

Given that the current Wh searches do not exclude the parameter space of interest, we are left 
with the VF + VF _ + ^ x final state produced from the decays of charged uj fermions with lifetime 
> 1 mm (see Fig. 8(a)). This is regarded as prompt decays by ATLAS and CMS, both of which 
require the transverse impact parameter to be |do| < 1 mm. Such prompt WW final states have 
already been discussed for Region (in). 

Region (vi): Displaced Wh + Displaced WW + $ T (diagrams in Fig. 8) 

This region displays a rich assortment of exotic signals, since both the charged and neutral uj decay 
through displaced vertices to xo, as in Fig. 8. The lifetimes of the charged uj states are in the 
range 1-50 mm, so they still decay before reaching the inner detector, rendering them safe from the 
searches for long-lived charged particles, as discussed for Region (i). On the other hand, the neutral 

6 A fraction of cjo will decay inside the inner detector, for which we refer the reader to section on Decays inside 
the inner detector or before in Region (vi). 
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(a) 



(b) 



Figure 9: Diagrams show 6-tagging using secondary vertex (SV) in processes where (a) the 6-quark 
(shown in red) is emitted from primary vertex (PV), and (b) the £>-hadron (shown in red) is emitted 
from the decay of a long-lived particle such as uo (shown in green). In both cases, the jet axis is 
represented by dashed black arrow, and the decay length vector L by the dashed blue arrow. Note 
that the decay length in diagram (a) coincides with the £>-hadron track. In both diagrams (a) and 
(b), do is the impact parameter for a particular track. Diagram (a) has positive signed impact 
parameter, while (b) has a negative signed impact parameter (refer to text for details). 

state, cjo, can have lifetimes as long as a few meters and can thus decay to Higgs and missing energy 
anywhere inside the detector. The relevant final states are displaced WW + $ T and Wh + $ T , 
where the W is displaced by at most ~ 5 cm, while the h can appear anywhere in the middle of the 
detector. The fraction of uo decays in various parts of the detector are plotted in Fig. 3(b). 

As discussed for Region (ii), the displaced W are difficult to search for due to the QCD back- 
ground in the hadronic decay channel and the inability to reconstruct the secondary vertex in the 
leptonic decay channel. Nevertheless, the leptons from displaced W decay can be used for triggering 
purposes, although trigger efficiency is expected to be low since the leptons may not be high-pr for 
our benchmark point. On the other hand, for the displaced Wh + $ T signal, the Higgs bosons can 
be used to reconstruct the corresponding secondary vertex in many cases. Below we discuss the 
collider signatures and relevant searches depending on where the ujo decays inside the detector: 

• Decays inside the inner detector or before: This section covers the parameter space with 
ujo lifetimes < 1 m, including the range 1-10 mm of Region (vi). For A3 between w 6 x 10 5 TeV 
and « 3 x 10 6 TeV, uo will decay to h + $ T predominantly inside the inner detector or the 
beam pipe. The Higgs boson thus produced will then decay to a number of final states: 

(a) h bb : This is the dominant decay mode of the Higgs boson, with BR(h — >► 66) « 58% 
for nth ~ 125 GeV. CDF and DO conducted searches [37] for long-lived neutral particles 
decaying to bb with lifetimes of ^ 2 cm and ~ 2.5-10 cm, respectively, in the context of 
hidden valley scenarios [4] , where the Higgs decays to a pair of long-lived neutral particles 
7r v , which in turn decay to 66, thus giving rise to h —> 7r v 7r v —> 2626. The applicability of 
these searches to our case is far from clear, however, because our h is not highly boosted 
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like 7r v , which should affect the efficiencies. Nonetheless, given that triggers already exist 
for highly displaced b jets (through high-px niuons produced from £>-hadron decays), it 
would be worthwhile to look into our displaced Higgs signals in the Tevatron data. 

(b) h n jets + X (n > 2): Two or more jets in the final states can arise from a number of 
Higgs decay channels, e.g., h — > gg, h -> 66, h -> WW* -> 4j or £v2], h — > ZZ* — > 4j or 
2j + X. If the cjQ decays to jets (through Higgs) between ~ 40 cm-100 cm of the ATLAS 
inner detector, i.e., outside the pixel detector but inside the silicon and transition radiation 
trackers, then the decays would appear as trackless jets. While currently no searches exist 
for this signature, it is our understanding that ATLAS has implemented [5] a signature- 
driven trigger for trackless jets in the hidden valley context. Such a trigger should also 
help selecting events in our model. 

(c) h £ + n jets + X: Single-lepton triggers may be used to select this channel, where the 
lepton comes from uj -> (h -> WW*/ZZ*) + $ T or u±/u>± -> W ± + $ T . In particular, 
the final state h WW* —> /iz/2j is similar to signatures of R-parity violating supersym- 
metric models where the long-lived neutralino decays to muons and jets with a displaced 
vertex. Hence, the ATLAS searches for such neutralinos [36], where the trigger requires a 
muon to have px > 50 GeV, should apply to our case with a similar sensitivity. 

(d) h —> ni + X (n > 2): Single-lepton or double-lepton triggers can be used to select such 
events. The multi-lepton final states can arise from ujo —> (h — > WW* / ZZ*) + $ T . A 
search for long-lived particles decaying to ££ + $ T was conducted at DO [38, 42], and 
we have checked that our benchmark point = 140 GeV is safe from this search due to 
weak exclusions at the Tevatron. The CMS experiment has searched for neutral long-lived 
resonances decaying to dileptons [39]. In our model, such displaced dilepton resonance 
is provided by the Z from uo — » (h ZZ*) + Our model, however, is safe from 
this CMS search because of a decreased sensitivity near the Z resonance as well as a 
large lepton px cuts used in the search. Searches for displaced and highly boosted muon 
jets [41] are not directly relevant to our scenario, as our muons will not be very boosted. 

(e) h 77 or ^yZ: Another smoking gun signature of this region of parameter space is the 
displaced 77 or 7Z resonance that reconstruct to the Higgs mass. Although there are 
existing searches for displaced single or two photons [35, 42], none of them seek a resonance 
from a diphoton pair, as they are primarily motivated by gauge mediated supersymmetric 
models. 

• Decays inside the calorimeter: There are currently no searches for long-lived particles 
decaying inside the calorimeter, except for stopped particles searches [43], which are not 
directly relevant for our scenario. ATLAS, however, has implemented triggers for long-lived 
neutral particles decaying inside the HC AL [5] . These triggers are motivated by hidden valley 
models but can be just as sensitive, if not more, to our displaced Higgs. Events are triggered 
by the following two signatures designed to characterize decays inside HCAL: (a) The decay 
products are confined to a small region inside the calorimeter, and (b) The ratio of energy 
deposition inside the ECAL to that in HCAL is very small. Depending on the luminosity, the 
benchmark point of our model can be sensitive to this search channel for A3 ~ 10 7 TeV and 
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Figure 10: Feynman diagrams leading to hh + final states in Region (vii) of the parameter 
space, where the charged u fermions predominantly decay to cjq. Particles too soft to be detected 
at the LHC are shown in orange. 

A 2 > 100 TeV. 

• Decays inside the muon spectrometer: Fig. 3(b) shows that a significant number of 
decays will occur in the muon spectrometer in this region of the parameter space. ATLAS has 
carried out a search [40] for hidden valley scenarios in which the long-lived neutral particles 
7r v decay inside the muon tracker. There is, however, an important difference between the tt v 
and our displaced particle, h. Namely, the tt v masses considered by the ATLAS search are 
~ 20-40 GeV, while the h mass is r^i 125 GeV. This implies that our displaced h will be far 
less boosted compared to the 7r v . Consequently, the trigger efficiencies would be lower for our 
case, as some fraction of the hits in the muon spectrometer due to our h would be associated 
to a wrong bunch crossing. A realistic estimation of the trigger efficiencies as well as the 
sensitivity of the ATLAS search for our scenario would require a detailed detector simulation, 
which is beyond the scope of a theory paper. 



Here, the charged uo states decay predominantly to ujo instead of xo with relatively short lifetimes 
< 10 cm. The lifetime of cjo, on the other hand, can be as long as a few meters. As shown in Fig. 10, 
the relevant final state is hh + $ T , where both Higgs bosons are macroscopically displaced from 
the primary vertex. The phenomenology of displaced Higgs has already been discussed in great 
detail for Regions (v) and (vi), although it should be pointed out that the efficiency for detecting a 
displaced Higgs signal in this region goes up by a factor of two, simply because there are two Higgs 
bosons. 



Region (vii): Displaced hh J r $ T (diagrams in Fig. 10) 



Region (viii): Pure $ T 
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In this region, the relevant diagrams are given by Fig. 10, but the neutral state cjq is so long-lived 
(> 10 m) that it decays predominantly outside the detector (see Fig. 3(b)). Similar to dark matter 
searches using $ T -\-X (X = j, 7, W, Z), this region can be probed at the LHC using jets, photons or 
W/Z bosons from the initial state partons. The relevant searches are : j [45], ^ T + 7 [46] and 

fi^ + W/Z [47]. For the J^x+j case, we generated the cjcj + Oj and uouo + lj (p J T > 80GeV) samples 
in MadGraph and showered/hadronized the partonic events in Pythia. The partonic cross-sections 
for the 0-jet and 1-jet samples were found to be 5.50 pb and 0.22 pb respectively at 7TeV. For the 
CMS, the efficiencies for passing $ T > 350 GeV, jet pt > 110 GeV and other cuts was found to be 
~ 0.2% and ~ 2% for 0-jet and 1-jet samples respectively. Even without performing any matching, 
we obtain at most ~ 70 events with a luminosity of 5 fb _1 whereas the 95% CL limit on non-SM like 
events placed by CMS from this search is 158 events. A similar analysis for the ATLAS experiment 
and for other final states with photons/ W/Z bosons yields the same conclusion that this region is 
not excluded by $ T + X searches. However, improved sensitivity is expected with 8 TeV data. 
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